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Abstract 

Using explicit perturbations of isotropic cosmological models which 
describe simple gravitational waves, an isotropic tensor having the al- 
gebraic symmetries of the Bel-Robinson tensor is derived as a model of 
cosmic background gravitational radiation and this is used to provide 
an answer to the question: in what sense can an energy-momentum- 
stress tensor similar to that describing the cosmic microwave back- 
ground radiation (neglecting anisotropics) be associated with an isotropic 
background of gravitational radiation? 
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1 Introduction 



It is well known (see, for example, [T]) that, neglecting small anisotropies, 
the cosmic microwave background radiation (CMB) in the universe is well 
described by a homogeneous perfect fluid which shares its 4-velocity with 
that of the isotropic matter distribution and has an equation of state giving 
the isotropic pressure equal to one third of the proper energy density. This 
paper is concerned with providing an answer to the question: in what sense 
is an isotropic background of gravitational radiation described by a similar 
perfect fluid to that of the CMB? To answer this question we make use of 

(a) simple explicit gravitational waves described as small perturbations of 
Friedmann-Lemaitre-Robertson-Walker (FLRW) cosmological models and 

(b) the Bel-Robinson tensor [2]. In section 2 we introduce notation and de- 
scribe a derivation of the perfect fluid energy-momentum-stress tensor for 
the CMB in a way which motivates our approach. In section 3 the corre- 
sponding calculation is performed on the Bel-Robinson tensor and in section 
4 a derivation of a perfect fluid energy-momentum-stress tensor for a back- 
ground of gravitational radiation is given. 



2 Cosmic Background Electromagnetic Radi- 
ation 

The metric tensor components (gab) corresponding to the FLRW cosmological 
models are given via the well-known Robert son- Walker line-element 

ds'^ = gabdx" dx^ = R^{t) f-^{dx'^ + dy^ + dz^} - dt^ , (2.1) 

with R{t) the scale factor and 

f = l + \{x' + y' + z')^ (2.2) 

where = 0, ±1 is the Gaussian curvature of the homogeneous hypersurfaces 
t = constant. The matter content of these universe models is a perfect fluid 
with the integral curves of the vector field d/dt = u'^d/dx"' the world lines 
of the fluid particles. Here x° = {x,y,z,t) for a = 1,2,3,4 respectively. 
The coordinate t is proper time and x, y, z are each constant along a fluid 
particle world line. The vector field with components u"' satisfies gabU°'u^ = 
u^'Ua = — 1 and is the 4-velocity of a fluid particle. The matter world lines 
are time-like geodesies, thus 

if = u^.bU^ = , (2.3) 



2 



with the semi-colon indicating covariant differentiation with respect to the 
Riemannian connection associated with the metric tensor Qab- Throughout 
this paper we follow the definitions and sign conventions of |3j and in partic- 
ular a dot on a tensor will indicate covariant differentiation in the direction 
of m". These geodesic world lines have vanishing shear 

O-ab = U{a;b) + U(aUb) " ^0 Kb , (2.4) 

with the round brackets denoting symmetrisation, hat = Qab + Ua Ub is the 
projection tensor and the scalar 6 denotes the expansion of the world lines 
given by 

e = ur^, = 3§. (2.5) 

The dot on R{t) denotes differentiation with respect to t. In addition the 
vorticity tensor 

(^ab = U[a-b] + U[aUb] , (2.6) 

vanishes in the case above. The square brackets here denote skew symmetri- 
sation. The energy-momentum-stress tensor of the matter is isotropic and 
thus has in general the perfect fluid form 

T'^'' = + (2.7) 

with proper density /x and isotropic pressure p functions of t only (and thus 
satisfying = = h^aP,b)- Einstein's field equations 

Rab = Tab — -^Tgab , (2.8) 

where Rab is the Ricci tensor obtained from the components of the Riemann 
curvature tensor according to the formula Rab = R'^acb, T = T"'a and a 
factor involving the gravitational constant has been absorbed into Tab for 
convenience, yield in the above case 

2R k 



We see from these that 



fi + P = ^{R^-RR + k) , (2.11) 
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and so it easily follows from f l2.9p and fl2.10p that 

A = -(/i + p)^, (2.12) 

e = -le'-^{^ + 3p) . (2.13) 

The first of these is the equation of conservation of energy along a fluid world 
line and the second is Raychaudhuri's equation. 

Electromagnetic test fields on the space-times described above have elec- 
tric and magnetic parts defined respectively by 

Ea = Fabu' and = *Fabu' , (2.14) 

where Fab = —Fba is a tensor field on the FLRW space-time having dual 
*Fab = IVab^'Frs with r]abcd = ^abcdy/^, 9 = det^Qab) and eabcd the four 
dimensional Levi-Civita permutation symbol. Maxwell's source-free field 
equations on the FLRW space-times read [3] 

^";a = 

and 

E"" + -6 E" 
3 

3 

To solve these equations we introduce a 4-potential a"" satisfying [1] 

a'^M, = 0, a% = 0, (2.18) 
from which the Maxwell tensor Fab is calculated using 

Fab = (Tb-a - Cra;b ■ (2.19) 

Now (12.151) and (12.171) are automatically satisfied while fl2.16p yields the wave 
equation (see |3]): 

a'^-"^,,= ^-{f^-p)a'^ . (2.20) 

Some simple sinusoidal, monochromatic electromagnetic waves on a A; = 
FLRW space-time satisfying (KT^ and iK2li\\ and thus f l2J[5|) -f l2T71) are 
constructed as follows (see [5] for families of test electromagnetic radiation 
fields on isotropic cosmologies): Let 

(f) = 6a(3n''x^ -T , (2.21) 



H^^a = , (2.15) 

= -r]^''^UbH,,d, (2.16) 
= r^^'^'ubE,,,. (2.17) 
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where = (x, y, z) with Greek indices taking values 1, 2, 3. Here dap is the 
three dimensional Kronecker delta and we note from fl2.ip that = R~'^gaj3. 
Also are constants and T is a function of t only satisfying 

n-n = 5„;3n"n^ = l and ^ = ' (2-22) 

respectively. With the metric tensor given via (12. ip with k = it follows 
that 

rt,a0,6 = O , (2.23) 

and so the hypersurfaces (/)(a;",t) = constant are null hypersurfaces in the 
FLRW space-time with k = 0. They will be the histories of the wave fronts 
of the electromagnetic radiation. As candidate solutions we take a potential 
1-form 

aa dx" = X 5a(3iB'^ COS + sin 0) rfx" . (2.24) 

Here B = (-B") and C = (C") are constant 3-vectors related to the 3-vector 
n = (ra") by 

B = b X n and C = b - (b • n) n , (2.25) 

for any 3-vector b such that b ■ b = 1 and for any real constant A. We note 
that 

n-B = n-C = B-C = and B-B = C-C = l-(b-n)2, (2.26) 

and we have the useful relation 

S^J.p {n" + (B ■ B)-i(fi" + C^)} = 6^, . (2.27) 

To verify that (12.241) satisfies (12.181) and (I2.20p we note that each of the two 
terms on the right hand side of (12.241) has the general form 

aadx" = SaF{^)dx'' , (2.28) 

for any function F. Then the equations in (I2.18P imply [1] 

s"Ma = 0, s^a = and s''0,a = O, (2.29) 

while (I2.20p provides us with the equations: 

rt,a0,6 = 0, (2.30) 

s'^''0,6 + ^0,/,s'^ = 0, (2.31) 
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s'^-'',, = ^(/x-p).V (2.32) 



Now each choice of Sa leading to fl2.24p is of the form 

Sa = {Sai3 a'^, 0) with a ■ n = , (2.33) 

where a is a constant 3-vector and (p is given by fl2.21l) . Thus with the metric 
tensor given via (12. ip with k = it easily follows that (12.290 and (I2.30p are 
satisfied. To verify that (12.311) and (I2.32p are satisfied requires 

0,/ = ^ , (2.34) 

p and /X given by (12. 9p and (I2.10p and the non-vanishing components of the 
Riemannian connection associated with the metric Qab with = 0, namely, 
r^4 = R/R6ai3 and F^^ = RRdap- We note that the propagation direction 
in the FLRW space-time with A; = of the histories of the waves is and 
we see from (12.340 that the null geodesic integral curves of this covariant 
vector field have non-vanishing expansion on account of the expansion of the 
universe and so the waves are not plane waves. The unit 3-vector n gives the 
direction of the wave propagation at any point in the 3-space t = constant. 
Substituting fl2:2lD into fl2:T9|) and then calculating (121^]) we find that Ea = 
{Ea, 0) and Ha = {H^, 0) with 

Ea = 4^«/3(-5'^sin0 + C^cos0) , (2.35) 
R 

Ha = -45a/3(5^cos0 + C"^sin0) , (2.36) 
R 

and in addition we have 

= —F4a = Ea and Fal^ = —Fpa = R^ajB-y H^ , (2.37) 

with eap-y the three dimensional Levi-Civita permutation symbol. The elec- 
tromagnetic energy tensor associated with these electromagnetic waves is 
calculated from 

Mab = \{Fac Fb" + *Fac * Ff,") . (2.38) 

Here Mab = Mba and = 0. Using ([23S])-(E3ZD and making use of the 
relation (I2.27p (in particular to simplify M^^ below) we arrive at 

M^, = ^(B-B)5^,5,^n"n^ (2.39) 

M^4 = -^(B-B)5^«n^ (2.40) 

M44 = ^(B-B). (2.41) 
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These components can be summarized economically as 



M,, = ^(B-B)0,,0,, , (2.42) 

with B ■ B given in terms of b and n by (12.261) . We now have explicitly 
how Mab depends upon the directions of b and n. To obtain a model of 
the background electromagnetic radiation in the k = FLRW universe we 
sum the electromagnetic energy-momentum tensors for all such wave systems 
labelled by the directions of b and n. Summing (i.e. integrating) Mat over 
the directions of b we find that 

M^, = ^S^aSufsW'nf' , (2.43) 

M^4 = -J^^f^-^"^ (2-44) 
M44 = (2.45) 

Now summing (i.e. integrating) over the directions of n we obtain the energy- 
momentum-stress tensor of the electromagnetic background radiation Aiat 
with components 

M,, = ^5,,, (2.46) 
= 0, (2.47) 

A^44 = ^, (2.48) 
where Cq = 327r^A^/9. This can be written in the form 

Mab = fJ'eUa Ub + Pehab , (2.49) 

with ^ 

= ^/Xe = , (2.50) 

and the subscript is used to reflect the electromagnetic origin of these quan- 
tities. Hence the model background electromagnetic radiation field that we 
have constructed in this way is a perfect fluid with fluid 4-velocity and 
equation of state fl2.50p . We note that A^"% = 0. 

We now generalize the construction above to include the FLRW universes 
with k = ±1. The end result is expected to be again the isotropic energy- 
momentum-stress tensor (I2.49p . Hence we can simplify the derivation by 
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calculating this tensor along any one of the integral curves of the vector field 
u"" and then extending the result to all curves on the basis of isotropy. For 
simplicity we work on the curve r = with = We will establish 

that given by fl2.33p . with given by fl2.2ip . in the neighbourhood of r = 
yields a solution of the equations fl2.18p and fl2.20p along r = 0. This is quite 
a remarkable fact which we will see later is not replicated in the gravitational 
radiation case. With it established we then have that fl2.49p will hold along 
r = and hence by isotropy will hold everywhere in the FLRW space-time 
with k = ±1. We note that (j) in fl2.2ip satisfies fl2.23p along r = with gab 
given via (12. ip and (12. 2 p with k = ±1. With / given by (12. 2p we have along 
r = 0: 

We shall require the components of the Riemannian connection associated 
with the metric tensor gab given via (12. ip and (12. 2p with k = ±1. These are 
listed in the appendix. Since now s° = (/^a"/i?^,0) and Sa = {6ai3a'^,0) we 
find that 



2f df f 
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From this it follows that along r = we have s"'.a = and {s"" ^a) ,bU^ = 
but {s°'-a),\ 7^ 0. Using (I2.5ip and the formulas in the appendix we find that 
along r = 

_ 2 d^f J_ a _ k_ _3k_ ^ 

-'""''^ ~ B^d^^ + ;^r„^,A^" - -^'^ - - -^SA . (2.53) 

Now with (s";J,A ^ along r = the equations ([222]), fl^^ and flOT]) 
hold (and are satisfied) along r = 0. However the wave equation (2.32) is 
modified to read 

Sa;d"' = l{l-i-p)Sa + K{s'',,),b, (2.54) 

along r = 0. This is trivially satisfied when a = 4 and thus we must have 

SA;/= ^(/i-p)sA + (s'^;d),A , (2.55) 

along r = 0. Direct calculation starting with s\ = and using the connec- 
tion components given in the appendix reveals that along r = 

^A;.' = + ^ + ^ j = -(/. - P)SA - ^ .A , (2.56) 

with the final equality arising from (12. 9 p and (I2.10p . Now (I2.55P is a conse- 
quence of (I2.53P and (12. 56 p . Since the passage from Sa (and thus from aa) 
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to Fab involves only first derivatives of Sa we obtain for Fab, and thus for 
Mab, along r = for k = ±1 the same expressions as those given by fl2.35p - 
fl2.37p and (12.421) in the case A; = 0. Hence summing the electromagnetic 
energy tensors as before with respect to the directions of the 3-vectors b and 
n yields the isotropic energy-momentum-stress tensor f l2.49p along r = 0. 
This construction can be carried out along any integral curve of the vector 
field u"" and thus the perfect fluid (12.491) with (l2.5Up models the background 
electromagnetic radiation in all cases k = 0, ±1. 



3 Cosmic Background Gravitational Radia- 
tion 

We consider now perturbations of the FLRW space-times which describe 
gravitational waves |6] using the gauge-invariant and covariant perturbation 
theory of Ellis and Bruni [7]. The perturbations we are interested in here 
describe pure gravitational waves [8J and are characterized by the existence 
of perturbations in the shear a""^ of the matter world lines and in the matter 
distribution itself in the form of anisotropic stress vr"^. The latter is a first 
order perturbation addition to the perfect fluid energy-momentum-stress 
tensor (12.70 satisfying vr"'' = tt^", 7i"-^Ub = = 7r"a. These perturbations 
generate perturbations in the "electric" and "magnetic" parts of the Weyl 
conformal curvature tensor Cabcd given respectively by [5] 

Eab = CapbgU^ U'^ aud Hab = *CapbqU^u'' , (3.1) 

where *Capbq = ^Vap^'^Crsbq- We note that this left dual of the Weyl tensor 
is equal to its right dual C*p^g = \rjbq^^Caprs- The Ricci identities provide us 
with [8j 

C^"';6 = , (3.2) 
together with the relationships between Eab, Hab and a"^, tt'^'': 

2 1 

Eab = -(^ab - -0 (Tab + -T^ab , (3.3) 

and 

Hab = -(Tia^''%)f9cU^ ■ (3.4) 

The following equations are obtained from the Bianchi identities [8]: 

E'^'.b = -^Ti'^'.b , H'^'m = , (3.5) 
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j^ah ^QE^b^ UrH^\.aV^>"^ = --(/^ + pK^ - -^"^ - -On""^ , (3.6) 
H- + BH-^u..EK,,^>"' = (3.7) 

Some simple sinusoidal, monochromatic gravitational waves on a A; = 
FLRW space-time are given by (see [8] for families of pure gravitational 
radiation fields on isotropic cosmologies) 

aahdx^dx^ = X R{t) 6af,Sp^ {{B^B'^-C^'C'') cos (P+iB^'C+CB") sin0}dx"dx'^ , 

(3.8) 

and 

R 

T^ah = '^-^CFab , (3.9) 

with A, 5", C° and as in ^TM . To verify that ([S3D and ([S3]) satisfy 

(13.71) we note that each of the two terms on the right hand sides of (13.81) and 

(13.91) has the general form 

(Tab = SabF{(i)) and Tlab = ^abF{(i)) , (3.10) 

with F arbitrary and Sab and Ilaf, symmetric, trace-free and orthogonal to u'^ 
on both indices. Following from (I3.2l) - (l3.8p the equations that Sab, ^ab must 
satisfy are [H| 

= , = , (3.11) 

s'^'A' + {]^<t>''-4 - s""' = , (3.12) 

and 

ab-4 _ Ig^ab + 1 1^ _ ^g2\ ab ^ _jja6 _ 2^j^^^ ^ 

'"'3 V 3^^ 3 9 / 3 ^ ^ 

For future reference it is helpful to note that Eab and Hab in (13. 4p substituted 
into (13. 6p results in a wave equation identical to (13.131) but with s"* and 11"'' 
replaced by a""^ and vr"'' respectively. Then substitution of fl3.10p into this 
wave equation, and utilizing the arbitrariness of the function F, yields in 
particular (I3.12p and (13.131) . Now writing 

Sabdx°'dx^ = R{t)aa(sdx"'dx'^ and Ilabdx°'dx^ = 2Raapdx'^dx^ , 

(3.14) 

with aa/3 = ttfia constants satisfying aajsn^ = = a^/j 5"^ with introduced 
in (I2.2ip it is straightforward to verify that (I3.1ip - (l3.13p are satisfied. 
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With aab and nab given by fl3.8l) and fl3.9p we find using f l3.3p and (13 ■4p 
that the only non-vanishing components of Eab and Hab are 



Eap = X6a^6p^{{C^'C'' - B^B") sin0 + {B^'C + C^'B") cos0}, (3.15) 
H^^ = X6^f,6^^{-{B^'C'' + C^'B'')sm(j)+{C^'C''-B^'B'')cos(j)} . (3.16) 

Alternatively on account of (13. ip the components of the perturbed Weyl 
tensor are now given by 

Cal3-y6 = — ^a/BX^^ScrExa , Cal3f4: = R^ajBuHa-i , C*q4/34 = -E'a/3 • (3-17) 

The analogue of the electromagnetic energy tensor of the previous section is 
the Bel-Robinson tensor [2] 

Mabcd = ■^{Ca^b'^ Ccpdq + *Ca^b'^ *Ccpdq) ■ (3.18) 

This satisfies M(^abcd) = Mabcd and M'^acd = 0. Using (I3.15l) - (l3.17p along with 
(I2.26P and (12.271) we find that Mabcd can be simplified to read 

Mabcd = 2\\B ■ B)20>_,0_,0,^ , (3.19) 

with B • B given by (I2.26p . Summing (i.e. integrating) this over all possible 
directions of b results in 

Mabcd = , r 4>,a4>,b(p,c(p,d ■ (3.20) 

15 

Finally summing over the possible directions of n we arrive at an isotropic 
tensor Aiabcd associated with the gravitational radiation background given 
by 

C2 1 

Mabcd = -^{UaUbUcUd + 2h^abUcUd) + -h^abhcd)} , (3.21) 

with cj = 2567r2AVl5. 

To extend this result to the FLRW space-times with = ±1 we follow 
the pattern of the electromagnetic case described in the previous section. It 
will be sufficient to establish f l3.2ip along the world line r = 0. By analogy 
with the electromagnetic case we might expect that the solution (13. Sp and 
(13. 9p might hold along r = in the case k = ±1 but this is not the case. The 
anisotropic stress perturbation (13. 9p can be retained but the matter shear 
perturbation has to be adjusted by the inclusion of a factor /^^^^ where / is 
given by (12. 2 p . The simplest way to see this is in terms of the variables (I3.14p . 
With or without the factor f'"^^^ multiplying Sab the field equations will be 
satisfied along r = except the final wave equation. This is because although 
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s -b = along r = and the derivative of the left hand side here along r = 
(i.e. in the direction u") also vanishes the spatial derivative h'j{s"'^ ■}j)-d does 
not vanish along r = 0, analogous to fl2.53p . The wave equation evaluated 
along r = now instead of fl3.13p reads 

_ 2 ( \ \q _ 1^2\ ab _ lu'^ia yb) ^ _f^ab _ ^Q^ab _ 

' 3 V 3^ 3 9 / 2^ 3 

(3.22) 

If we put Sap = /" R(t) Qap (with a view to determining the real exponent 
n) and IIq,^ = 2R a^p then along r = we find that 

= ^^^^^^"/^ ' (3-23) 

and fl3.22p reads explicitly 

(R 3i?2 2k\ 3fc , , , 

-[-^ + ^ + ^2)''^^- ^("^ - = ' (3-24) 

along r = 0. Direct calculation of the first term here, using the formulas in 
the appendix, yields 



s 



al3;d 



3nk R 3R^ 3k 



2R'^ \R i?2 



along r = 0. Now fl3.25p in fl3.24p results in the requirement that n = 
—7/3. Since the perturbed Weyl tensor described by Eat, Hab given by (13. 3p 
and (13. 4p only involves first derivatives of aah.T^ab the perturbed Weyl tensor 
calculated along r = will in the case k = ±1 be given again by (I3.15P and 
(I3.16p . leading once again to (I3.2ip along r = and thus, by isotropy, along 
any integral curve of the vector field m". 



4 Energy— Momentum— Stress Tensor for Grav- 
itational Background 

The energy-momentum-stress tensor Aiab in (I2.49P describes macroscopi- 
cally the CMB radiation. It has been derived starting from an electromag- 
netic energy tensor Mab in (12.381) . J^ab has the same algebraic symmetries as 
Mab in (I2.38P (symmetric and trace-free) but Aiab is 'fiot an electromagnetic 
energy tensor. If it were, and the corresponding electric and magnetic fields 
Ea and Ha were given in general by (12.141) . then they would have to satisfy 
the algebraic equations 

EaHb = EbHa , (4.1) 
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and 



EaE, + HaH,= -{E,E'^ + H,H') h. 



(4.2) 



because Aiab includes no energy flux relative to u'* and no anisotropic stress 
(see [3], p. 14). These equations imply that Ea = = Ha- In similar fashion 
the tensor Mated in fl3.2ip has the same symmetries as the Bel-Robinson 
tensor flS.lSp (symmetric and trace-free on any pair of indices) but it is not 
a Bel-Robinson tensor. If it were a Bel-Robinson tensor then it would be 
derivable from electric and magnetic parts Eat and Hat of a Weyl tensor 
given in general by fl3.ll) . These parts would then have to satisfy analogous 
equations to fl4.ll) and fl4.2p namely, 



Eab Ecd + Hah Hcd ^ ^ (2 ^ac hd + 3 had he " 2 hab Kd) , (4.4) 

with Q = EabE"'' + HabH"'' (see [9], eqns.(5.1) and (5.2)). But these equa- 
tions imply that Eab = = Hat- 
In this paper we construct isotropic (with respect to the congruence tan- 
gent to m") energy-momentum-stress tensors by summing over waves propa- 
gating in all directions. In this context it is perhaps helpful to draw attention 
to the fact that no similarly isotropic electromagnetic field, described by Ea 
and Ha, exists (since these vectors are orthogonal to u"") but tensors quadratic 
in these quantities such as EaEb have non-vanishing isotropic forms which 
are proportional to the projection tensor defined following fl2.4l) . A similar 
comment applies to the tensor fields Eab and Hab- 

Mated in (13.211) is a gauge-invariant tensor field on the FLRW space-times 
which is small of second order (since it is quadratic in first order perturbations 
of the FLRW space-times), has dimensions of (length)"^ and has divergence 
given by 



Hat Ecd + Hcb Ead — Had Edj + Hcd Eab 



(4.3) 



and 



Q 



M 



ated 




On the other hand we find that 



{UaUbM''"'^),d = , 



(4.6) 



where by fl3.21j) 



UaUbM 



abed 




{UcUd + -hcd) - 



(4.7) 
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This tensor is symmetric, trace-free and divergence-free. It has the attributes 
of an energy-momentum-stress tensor except for the fact that its dimensions 
are (length)"^. However using the scale factor R{t) we can define an energy- 
momentum-stress tensor with dimensions (length)"^ by multiplying (14. 7p by 
Rl to give 

= flgU'^u'' +Pgh^'' , (4.8) 

with 

P. - i., ^ H . (4.9) 

for some Rq = -R(to) 7^ and the subscripts on p and reflect the gravita- 
tional wave origin of these quantities. In this way we associate an energy- 
momentum- stress tensor of the CMB type with an isotropic background of 
gravitational radiation. It is a second order perturbation of a background 
isotropic energy-momentum-stress tensor such as that describing the CMB, 
neglecting anisotropics. It would be interesting to incorporate it into the 
EUis-Bruni perturbation theory at second order. The energy-momentum- 
stress tensor (14. 8 p will act as a source of second order perturbations of the 
isotropic cosmologies. Anisotropics in this cosmic background gravitational 
radiation will thus be third order. 
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A Useful Riemannian Connection Com- 
ponents 

The non-vanishing components of the Riemannian connection associ- 
ated with the metric tensor Qab given via the line-element (12.1 p with 
(1^ are 

_pi _ pi _ pi _ f-i^ pi __f-i^ pi -_f-i^ (\u 

■p2 _ 1^2 _ -pi _ J--1 
-i22-ill-i33-/ ^'^12--/ 5^'^23--/ (A.2j 

_p3 _ p3 _ p3 _ f'l^ p3 _ _„-l^ p3 __f-l^ (\o\ 
i33-ill-i22-/ 5^'^13- Po5^'^23- / ^^'l^'^) 

r^^ = /2i?/-2 5,^ , (A.4) 

and 

r^4 = |5a^- (A.5) 

From these we see that along r = the only non-vanishing components 
are F^^ = RR5aj3 and F^4 = RjRdap- A further useful consequence is 
that along r = 

a,. = -^5«/.- (A.6) 
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